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Levchenko and Solov’ev (1972, 1974) have developed a stability theory for 
space periodic flows, assuming that the Floquet theory is applicable to 
partial differential equations. In the present paper, this approach is 
extended to unsteady periodic flows. A complete unsteady formulation of 
the stability problem is obtained, and the stability characteristics over 
an oscillation period are determined from the solution of the problem. 
Calculations carried out for an oscillating incompressible boundary layer 
on a plate showed that the boundary layer flow may be regarded as a 
locally parallel flow. 
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STABILITY OF AN OSCILLATING BOUNDARY LAYER 
V. Ya. Levchenko and A. S. Solov'yev 

1. The problem of stability of nonstationary, and, in particu- /10 * 
lar, periodic flows has long attracted the attention of the research- 
ers because of its high theoretical and practical importance. Experi- 
mental observations of the transition to turbulence in some periodic 
flows revealed a number of interesting phenomena. In attempting 
to explain them, the researchers have turned to the theory of hydro- 
dynamic stability. Early studies were based on the concept of 
quasistationary flow, which studied the stability of momentary 
velocity profiles [1, 2]. The quasistationary approach is apparently 
applicable if the change in velocity profile occurs fairly slowly 
in time, i.e., when the characteristic time for change in the main 
flow is considerably lower than the period of natural oscillations 
in the flow. However, the quasistationary approach is not satis- 
factory for many real situations, and could even lead to erroneous 
results [3]. Experiments on the transition to turbulence in 
oscillating boundary layers show that nonstationariness has a sig- 
nificant influence on the transition when the frequencies of the 
applied oscillations are on the order of unstable natural frequencies 
of the stationary boundary layer (Tollmin-Shlip.hting waves) [4-7]. 

There is a problem in applying the quasistationary model to these 
conditions . 

Publications [8,9] have developed a nonlocal theory for 
stability of 'spatial-periodic flows which is based on a hypothesis 
of applicability of the Flock theory to p.-nM 1 a I <1 1 riVrcii I 1 ;i I 
equations. The correctness of this theory was confirmed 
by direct experiments [10]. In this work, the methods developed 

in [8,9] were extended to nonstationary periodic flows, i.e., an 
^ » 

Numbers in right margin indicate pagination in original text. 


1 


attempt was made to set up and solve the problem of stability of 
these flows in a "complete, nonstationary situation," thus defining 
the stability characteristics for the oscillation period. Specific 
calculations were made for an oscillating boundary layer of 
incompressible liquid on a flat plate; the flow in the boundary 
layer in this case was considered to be "locally parallel." 

A number of researchers in recent years, independently of each 
other, have abandoned the quasistationary approach and conducted 
studies on stability of oscillating flows in a "complete nonstation- 
ary formulation. Publication [3] solves the problem of stability 

of a layer of viscous liquid with free surface that is put into 
motion by the periodically oscillating lower boundary in relation 
to long-wave perturbations. Grosch and Salven [11] studied the 
stability of a Poiseuille modulated planar flow, von Kerczek 
and Davis [12] studied the stability of a finite Stokes layer 
between two parallel plates , one. of which was brought into harmonic 
oscillation in its plane. These publications are also based on the /II 
aforementioned hypothesis. The technique for computing the 
stability characteristics in these works varies and differs from 
that proposed by the authors. The authors do not know a solution 
to the problem of stability of an oscillating boundary layer in 
a "complete nonstationary formulation." 


2. The behavior of small two-dimensional perturbations in a 
flat flow of viscous incompressible liquid is described by flow 
function t|^(x, y, t), which satisfies the linearized vortex equation 



djM\ if) 
d (X, y) 


d(W. 1 ° ) 
d (x, y) 


■aa r , 



< 1 * 

du 2 * 


( 2 . 1 ) 
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the examined case is a function of time t and the coordinate y, 
perpendicular to the direction of the flow; R--Reynolds number 
which is constructed using characteristic dimensional quantities 
of length and velocity. The function ijj° does not depend on the 
longitudinal coordinate x, therefore, one can look for the solution 
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to equation (2.1) in the form 


V U, t) —e !axX V (t, y). 


( 2 . 2 ) 


We will examine the spatial growth in small perturbations, 
considering that a— a r -\-iai ' is a complex constant. In t his case, 

the flow is stable if ^>0, and unstable if 0. 

Assume that the main flow \|>°(t,y) is periodic for t with period 
2ir/e. Assuming that one can use the Flock theory for ordinary 
differential equations for the partial differential equation 
(2.1), we will look for the solution to (2.1) in the form 

M (2.3) 

where w--real constant, where spatial growth of the perturbations 
is examined, while $(t,y) --periodic function t with the same 
period 2 tt/e as the coefficients of equation (2.1). Some sub- 
stantiations for this hypothesis are presented in publications 
[3, 11, 12]. By substituting into (2.1) functions ip°(t,y) and 
$(t.y) in the form of infinite Fourier series for t with coefficients 
that depend on y, 

30 CO 

V ( t . y) = 2 K (y) e iM -I- S hi (y) c iM , (2.4) 

ll.--.--Q 

oo 

{/)= 2 q>« (it) c'"" 

— oo 

(2.5) 

and equating the terms with the same exponents, we obtain an 
infinite connected system of linear, ordinary differential equat Lons 

L„(p„ SS (u — c -f n j A„ — U"<p„ — ( rp 1 , v — 2a 2 (f’„ -f a\ n ) = (2.6) 

OO 

= 2 + /t/ { IU (p/H-ft — h k A„~k — /- 1 , 1 +fe) ( n = 0 , ± 1, . . .). 

k~ l 
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Here the following designations are introduced 

U = 2 Re h'o , ® = — ac < 4,< = <P* — a °" f Pn- 

The asterisk designates the complexly-conjugate quantity, the 
apostrophe designates the derivative for y. 

If the boundary conditions for <|> are uniform, then with assigned /I 
a) , R and e we have a problem for eigen . values, and the flow is 
stable or unstable depending on the sign of The i|>^(t,y) flow 

which is studied for stability is known before the problem is 
solved, therefore, in the sums for k in expressions (3.4) and (2.6) 
one can be limited with the necessary accuracy to a finite number 
of terms in the Fourier series. The main difficulty in solving 
system (2.6) is that system (2.6) is coupled. We assume that if 
series (2.5) is also broken, i.e., we ignore in (2.6) the terms 
containing <p, +1 , cp s+2) cp _ s _ 2 then, starting from a certain 

fairly large number n = s, where s is a whole positive number, 
the eigen values a found will not depend on s essentially with 
the necessary accuracy. The system of equations that is thus 
obtained has a finite order, and can be integrated. In particular, 
with k = 1 , it looks like 

!/i*Pn = fp„ — i -{- 4l (pn+l — Att+U X 11 — — S,...,s). (2 7 ) 

The plan for integrating systems of this type is presented in 
detail in [8]. It is clear that the effectiveness of the described 
method for the solution depends significantly on the number of 
terms s taken into consideration in series (2.5). The use of 
the method is justified, if the value s is small, i.e. , does not 
exceed two or three. The value s is selected while solving the 
problem and depends on the type of junction ^(t,y). However, 
calculations show that with a reduction in the frequency parameter 
e, starting from certain values, one has to take into consideration 
a greater number of terms in the Fourier series. This results in 
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an undesirable increase in the order of the corresponding equation 
system, increases the time for numerical integration of the system 
on a computer, and the described direct method for solving the 
problem becomes ineffective. 

On the other hand, when e is small, equation (2.1) contains 
coefficients which are slowly changing functions of time t, in 
that sense that their derivatives are proportional to the small 
parameter e. Asymptotic methods can be used in this case to 
solve equation (2.1). In this work, as in publication [9], the 
method of many scales is used in the process of constructing the 
asymptotic method to solve the problem in a complete nonstationary 
formulation. 

3. In order to simplify the calculations, we will examine a 
case frequently encountered in the applied cases where t|/^(t,y) 
looks like (in (2.4) k = 1) 

lP(t,y)=$o{y)+e°{li{y)e ,t, +h t {y)e~ M }, ( 3 . 1 ) 

where i|) 0 =2Re/io, 4i = e°4, . while e --a certain small amplitude para- 

" o 

meter. The parameters e and e satisfy the conditions 

e<l, e°<l. (3.2) 

We will introduce into (2.1) and (3.1) a new time scale T = Et, 

and, by using the formalism of the method of many scales, we will 

replace the derivative for t in (2.1) of the tyne A » 9 ■ - d 

dt^ dt~r R W 

In the resulting equation, the equations will not depend on x and 
* 

t, consequently, one can search for the solution in the form 

*'{x,y,T,t) = <Hy,T)eK”-', (3>3)> 

where the following equation is obtained for the function $(y, T) 
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(3.4) 


i<i 0 -f- e 


dT 




ia 


flip 0 ' 

dy 


)(■&-«) 




ia, 


dy* 


(I) = 


d* 

Oy- 


a 2 


After this equation has been supplemented with uniform boundary /13 
conditions, we obtain a problem for eigen values to determine a. 

The method for solving this type of problem is presented in detail 
in [9], therefore, we will present below only a brief description 
of the method and derive the main equations for the examined problem. 

The coefficients of equation (3.4) are periodic functions T 
with period 2 tt . We have a solution to equation (3.4) in the class 
of functions which are periodic for T with the same period 2 tt , 
and require that the two small parameters e ° and E are proportional 
to each other with coefficient of proportionality k, i.e., 


As indicated in [9] and discussed in detail in [13] in the example 
of a model equation, imposition of condition ( 3 . 5 ) is the only 
possibility of avoiding a quasistationary approximation . Decom- 
posing a and <j>(t,y) into asymptotic series for e 

'OC = CCo + e CX I — |— B 2 CC2 — ? 3 OC3 -f- . , . , 

<D(7\ r/) = ® 0 (T, y)+e'(D l (T, ij)+ £ 2® 2 (T, //) +e>(J) 3 (T, (/)+..., (3.6) 


where all <K periodic functions with period 2 t r , substituting 
(3.6) into (3.4) with regard for (3.5), and equating the terms 
with the same degrees of e, we obtain a series of partial 
di fToronti al oqnnf.1 onr. , I'i-oiii which i l.c ,v ,, j , ■ i , ,, r ,, n . 

(3.0), constants u i and iuncLion ^ are successively determined. 
The terms of the zero order yield the equation 


I<Do=0, 


(3.7) 

^ --operator of 


where L = ( - ico + i* 0 ^ 0 ) - a 5 ) - - -1 (JL _ 

Orr- Sommer field . One can search for the solution to this equation 



in the form 


®o {T, y)=<fo(y)B(T). 


(3.8) 


The eigen value. aQ and function <f>o(y) are found from the Orr 
-“Sommertield equation uniform boundary conditions for <)>q, while the 
B(T) function is determined from a nonuniform equation obtained 
by leveling the terms with e* 

L(D l = -kB(T)e‘ T F n (y)-kB(T)e-‘rp u {y)^ 

dB r, V ^ y / 

~ If (!/) ~ a i B ( T ) F 14 (y), 


where 


Fu(y) =itxo(h' A(fo—h'"<fo) , 

F n{y) == icLo(h t ' Acpo h cpo), 

F\z{y) = Acpo. Fn(!/)=*-2a 0 (- ico -f ia 0 %) - hp"'cp 0 +(i¥ 0 


(3.10) 


This equation, supplemented by uniform boundary conditions for 
<J >2 has a solution only when the condition of orthogonality is' 
fulfilled 


dB 

dT 


00 CO 

J cp+/V</ + ke!TB (T) | < + ke ~ iTB (T) I it f 1 Jy + M ( 7 1 ) X ( 3 . 11 ) 

F u dy = 0. 


The solutions to the equation, conjugate to the Orr- Sommer field /I A 

equation are designated here by <(>q. By using the requirement 
for periodicity of the B(T) function, it is easy to show [9] 
that = 0, while the B(T) function looks like 

B(T) =Bo(/')cxp[//f(FiiC? ,T — Vi 2 e~ ,T )], /o 19) 


where 


oo 



oo 

J vtFiM 


0 


(m = 1, 2), 
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while Bq ( k) --constant that is not defined in framework of the 
linear theory. We search for the solution to equation (3.9) in 
the form 


(3.13) 

®><0. r) <r) (£/) (^) ] +fli 

where the functions <J >2 and <f > ^ are found from the nonuniform 
Orr-Sommerf ield equations: 

icp m ~/j[ — F i„, (y) d-Timf i»(y)] 1, 2) (3.14) 


with uniform boundary conditions for <j> . In order to determine 
the function B^(T) it is necessary to examine the terms with 
The equation for ^ looks like 


L<S> i =*—kB l (T)e IT F 11 (y)—kB i (T)e~ iT F 12 (y) - d *LF u {y)-B (T) e" F 2l (y)~ 

-B(T) e" F 22 (y) —B(T) F a (y)-B ('/') e-* F u (y) - 
~B(T) e~ 2IT F 2 5 (//) —a 2 B (T) F u (y). 


(3.15) 


Here 

F 21 = k (fff - r n fi3 ), F n = iF\'l 

F n = k (fft - r ls fi , 3 ) .+ F?? - r 11 / r( i 2 3 ) ), " 

F 2i = - IF®, F n == k [F® - T n F?J ). 


(3.16) 


From the condition of orthogonality, B-^T) is determined, while 
the requirement for periodicity yields 


a 2 — 


j <po F ia‘iy 
0 

J <Po F udy 


(3.17) 


Continuing the described algorithm, we obtain a^, etc. quite 
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analogously. 


a 3 



no 


f <Po ~ F ud(! 
o 


(3.18) 


where 


F ni = k [F® - T n F S 3 3 > + F\ V - r n F\$). (3.19) 

We note that in order to obtain the functions , in (3.16) 

and (3.19), it is necessary to replace the function by the 
function <j>^ (m = 1 , . . . ,4) . The functions <t>-^ and o are found from 
(3.14), while <j>g and <j>^ are determined from the equations 

4cp 3 = - y ikT n (- F u + T 12 F U ) - (F Si + r 22 F 13 ), (3.20) 

L< fi = 1 ik r, 3 (- Fjj + r u F la ) - (f m + r 24 f 13 ) 

with the corresponding uniform boundary conditions. Here 

oo 

J { ft F 2n dl J 

r 2'> = ^ — ( /Z==1 5). . (3.21) 

j’ To'^13 dy 


Thus, by analyzing (3.18) - (3.20) with regard for (3.14) 
and (3.16) and (3.21), and taking (3.6) into consideration with 
accuracy to terms on the order of e^, we find 

a=ao4-e 2 ^ 2 Co+e 3 (4 2 Ci+/e 4 Q)+0(8 4 ). (3.22) 

Here C n (n = 0,...,2) are known constants which depend on a Q , R, w 
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and do not depend on e and k. 


4. The mode of flow in the oscillating boundary layer with 
rate of the external flow U = 1 + e°coset is characterized by 
the parameter £ = The dash notes dimensional quantities. 

Here x-- longitudinal coordinate, ^--average rate of external 
flow, ^--cyclic frequency of applied oscillations. With t, < 1, 
the flow is considered to be quasistationary, with c > 10 it is 
high-frequency; the values 1 < x, < 10 correspond to the intermediate 
mode. Ackerberg and Phillips obtained asymptotic solution to 
the equations of the oscillating boundary layer on a flat pi ate 
for e -> co which is correct to values x, > 4 [1.41. This solution 
was used in calculating the stability characteristics. 


Figure 1 illustrates the neutral curves, figure 2 shows 
the curves for the coefficients of _perturbation amplification 
with fixed frequency parameter F = =3- — 110 • 10~ 6 , , where <d- -dimensional 

frequency of perturbation ,v--viscosity , while the amplif icatiion 
coefficients are designated by ina = - 


a ( -§dR. 


Curves 1 were 


computed for a stationary boundary layer (Blasius profile) . Curves 
2 are the results of calculations for the amplitude of applied 
oscillations e° = 0.1, and dimensionless frequency f # 45- 10" 6 . 


The calculations were made by both aforementioned methods. 

In that region of frequencies where e is fairly small and it can 
be viewed as a small parameter of the problem, the asymptotic 
calculation method was used (we note that c — FqR diminishes with 
a decline in R and F Q ; R- -Reynolds number constructed using the 
characteristic dimension g __ ^ j _ The neutral curve 2 in 

figure 1 was calculated both by the direct method described in 
section 2 (two terms in the Fourier expansion were considered, 
s = 2), and by the asymptotic method, using formula (3.22) (the 
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Figure I . 



corresponding points of the neutral curve are designated by light 
circles). The results show, the good agreement of the presented 
methods . 

_ r 

The dimensionless frequency Fq = 45 x 10 was selected to 
be close to the el S en frequency of the stationary boundary 
layer observed near the point of transition to turbulence [15] . 

a 

When acoustic oscillations are applied (amplitudes of the 
oscillations are small) with frequencies that coincide with 
unstable eigen ' frequencies of the stationary boundary layer that 
correspond to the internal region of the neutral stability curve, 
a significant increase was observed in the experiments in the 
Reynolds transition number, i.e., destabilization of the boundary 
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layer [4-7]. Since in boundary layer type flows, the 

transition to turbulence is due to the instability in the laminar 
flow in relation to small perturbations, one should expect the 
destabilizing influence of oscillations in the external flow on the 
stability characteristics of the boundary layer. However, the 
calculation results were unexpected and seemingly contradictory 
to the known experimental facts. For the value of the amplitude 
of oscillation in external flow = 0.05, the calculations indi- 
cated that the influence of nonstationary flow on the characteristics 
of its stability is essentially missing. With an increase in the 
modulation amplitude of the external flow, the flow stabilizes 
in the boundary layer: the curve of neutral stability for = 0.1 

(curve 2 in fig. 1) is shifted into the region of large Reynolds 
numbers; the slowing down of growth in perturbation is visible in 
this case in figure 2 (curve 2). 

Direct experimental verification partially confirmed the 
correctness of. this theory: with low amplitudes of applied 

oscillations, no change was noted in the stability characteristics 
of the boundary layer ; destabilization in the oscillating boundary 
layer observed in the experiments is not explained by the increased 
flow stability, but by other reasons [16]. 

Stabilization of the boundary layer with "high-frequency" 
oscillations in the external flow with fairly -high amplitude 
has not yet been experimentally confirmed; there are no experi- 
mental data for these conditions. The possibility of stabilization 
has basically been excluded. Stabilization was experimentally 
observed for flow in the boundary layer under the influence of 
sound waves transverse to the direction of the flow. The stabi- 
lizing effect of modulation on the flow of other types is known 
[11]. The results of this work can be qualitatively compared with 
the results of von Kerczek and Davis [12]. With high-frequency 
oscillations, all the influence of nonstationariness on flow in 
the boundary layer is manifest in a narrow band near the surface, 
the so-called Stokes layer. The calculations of von Kerczek and 
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Davis indicated very high stability of the finite Stokes layer, 
and the stability improves with an increase in the oscillation 
amplitude . 

The experiments of Obremski and Fejer [17] studied the 
effect of low-frequency oscillations on the transition in the 
boundary layer. Destabilization was observed in the flow with 
values of parameter p eyvs _j£ ( larger than 2,6 x 10^. With Re^g < 

2.6 x 10^ , there is no influence of nonstationary flow on the 
transition in these experiments. In the calculations, whose 
results are discussed above, the Re^g values did not exceed 2500. 

In figure 1, the cross designates the neutral points that were 
calculated for Re^g = 10 , Fq = 3 x 10 , which corresponds 

roughly to mode No. 2 in the experiments [17] , in which significant 
reduction was observed in the Reynolds transition number. Unfor- 
tunately, we did not succeed in advancing into the region of smaller 
Reynolds numbers with the employed Velocity profiles or con- 
structing the entire neutral curve, since the values of the c 
parameter in this case became smaller than 4. However, there 
was a noticeable expansion in the range of unstable frequencies, 
which indicates destabilization of the flow. Of course, it is 
impossible to draw definite conclusions from two points, but it 
is possible that in contrast to the case of "high frequencies", 
destabilization of the flow at low oscillation frequencies of 
the external flow is explained by increased instability of the 
nonstationary velocity profile. Further studies of this question 
are needed. 
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